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Phase structure and chiral limit of compact lattice QED with Wilson 
fermions * ' 
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We study the phase structure and chiral limit of Ad compact lattice QED with Wilson fermions (both dynamical 
and quenched). We use the standard Wilson action (WA) and also the modified action (MA) with some lattice 
artifacts suppressed. We show that lattice artifacts influence the distributions of eigenvalues Xi of the fermionic 
matrix especially for small values of Xi . Our main conclusion is that the chiral limit of compact QED can be 
efficiently located using different techniques. 



1. Introduction 

The lattice formulation of QED is not unique. 
One has to decide on a physical ground which 
version of QED is realised in nature if different 
lattice versions of QED do not belong to the same 
universality class. The old physics, i.e., known 
from experiment, has to be reproduced. 

When we consider QED as arising from a sub- 
group of some non-abelian (e.g., grand unified) 
gauge theory we are necessarily dealing with the 
compact version. Our choice in this work is a 
compact formulation of QED. 

In the theory with Wilson's fermions chi- 
ral symmetry is broken explicitly, and, presum- 
ably, can be only restored by fine-tuning the pa- 
rameters in the continuum limit if we are dealing 
with a meaningful lattice discretization. What 
one can expect at nonzero spacing is that at some 
k c = k c {P) a so-called partial symmetry restora- 
tion takes place |2|,|| when the Wilson mass term 
and ordinary mass term cancel at zero momen- 
tum in certain vertex functions. If so one can 
approach continuum limit and chiral symmetry 
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restoration along the line n c ((3). Another ques- 
tion is if in the continuum limit of our theory the 
chiral symmetry is realised explicitly or is spon- 
taneously broken. 

2. Actions and order parameters 

The modified lattice action Sma(U, 4>, 4>) for Ad 
£7(1) gauge theory (QED) is 

Sma = P'S g (U) + Sf(U,$,iI>) + 8Sg(U). (1) 

In eq.(|l|) Sa(U) is the standard plaquette (Wil- 
son) action for the pure gauge U(l) theory, and 
the additional term SSq suppresses lattice arti- 
facts (i.e., monopoles and negative plaquettes). 
The fermionic part of the action Sf(U, ip, ip) is 

N f 4 
f=l x,y s,s' = l 

M = \ + k-M{U), (2) 

where M. is Wilson's fermionic matrix, Nf 
is the number of flavours and n is the hopping 
parameter. The first two terms in r.h.s of eq.(jl]) 
make up the standard Wilson action Swa- 

In our work we used both Swa and Sma with 
Nf — 2 for dynamical fermions. Apart from (ipip) 
and (V^75?A) we calculated the pion norm (II) jj] 

U ( U ) = jj ■ J2 ^(M-^M-^ B ) , (3) 



2 



where L is the lattice size. Its advantage is that 
it appears to be a very sensitive quantity in the 
'critical' region. Introducing eigenvectors g n (s, x) 
of with eigenvalues fi n : ^Mg n = M« - 5n , 

one can easily obtain a spectral representation 

n = FEiEi3«M)i 2 • w 

n ' s 

Following the common practice we identify here 
the chiral transition with the appearance of zero 
or near to zero eigenvalues of the fermionic matrix 
M.. Evidently, an eigenstate of M. with eigen- 
value zero is also an eigenstate of 75.M. So, the 
presence of configurations which belong to zero 
eigenvalues of M. gives rise to poles in II . 

3. Conjugate gradient (eg ) method 

To locate k c (P) one can use the convergence 
rate of the cg-method. This is the iterative 
method of solving a system of linear equations 
D ■ X = ip, where D is a hermitian n x n matrix 
and <p is an input vector (D = Ai'Ai in our case). 
The convergence of the cg-method should be con- 
trolled by the condition number £ = X max /X m i n , 
where X max and \ rn in are maximal and mini- 
mal eigenvalues of D. Close to k c the minimal 
eigenvalue of Ad^Ai is small and is supposed to 
be ~ (l — k/ k c ) ■ 

We observed that at large enough number of it- 
eration steps N cg (N cg > N ) the average residue 
(R) = (R){N C g) behaves as 

(R)=C- exp(-a • N cg ), a = In ^±1 (5) 

V? — 1 

independently of the distribution of eigenvalues 
Xi provided n is large enough. To check it we 
generated D with different (uniform, gaussian, 
double-peaked) distributions of eigenvalues and 
given Xmin and X max . The components of the 
initial vector Xq and of (p were chosen every 
time randomly with gaussian distributions. 

It follows from the above that, for the inversion 
of M^Ai , the (N^g 1 ) required for convergence to 
some small but fixed R will behave as (N^ 1 ) ~ 
1 — k/k c . In Fig.l we show the dependence of 
(iV^ 1 ) on k at (3 = 0.8 for quenched MA 
(qMA) for different L. 



Figure 1. {Nz^){k) for qMA for different L. 



At k ~ k c the data fit nicely to straight lines 
giving reasonable estimation of n c . The volume 
dependence becomes rather weak for larger L. A 
similar picture was obtained for the case of dy- 
namical fermions. 

4. Phase diagram and chiral limit 

For WA with dynamical fermions (dWA) at 
P < (3q ~ 1.0 thermal cycles with respect to k 
or p have a typical hysteresis behaviour for tj)tp 
and plaquette □. Time histories (TH's) of ^ip 
and plaquette □ for different starts show the ex- 
istence of metastable states. So, we conclude that 
for dWA there is a 1st order phase transition (PT) 
line from (/?;«) ~ (l.;0.) to {(3; k) = (0.;0.25) 
which is in agreement with 

After the suppression of lattice artifacts (i.e., 
for dMA) this line disappears (see also ||). 

At P > (3q and k < k c (P) the system with dWA 
is in the Coulomb phase. The photon correlator 
r(r) is well consistent with that corresponding to 
a zero photon mass. For /3 > Po and k > n c {P) 
the correlator T(t) shows a tachyonic-type be- 
haviour (mi < 0). Thus, we can conclude that 
there is a higher-K phase (or phases) differing 
from the Coulomb phase. 

For quenched WA (qWA) at P < p TH's of 
II (as well as of ipip and V'75'0) show very sharp 
peaks at n ~ Kc(P) which means the appearance 
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of small eigenvalues A^. It is worth noting that 
for WA at (3 < (3q those peaks do not disappear 
at k > k c but instead become even more strong. 

For WA at (3 > (3 and for MA at any (posi- 
tive) (3 the dependence of TH's on k changes dras- 
tically. We don't find peaks of comparable ampli- 
tude (~ 10 4 ) but, nevertheless, in some 'critical' 
region n <~ n c TH's of II become much more 
rough than at smaller or larger values of n. As 
far as for every configuration II is the arithmetic 
average of 4L 4 terms corresponding to 4L 4 dif- 
ferent eigenvalues one can conclude that rather 
small A, appeared. 

The (renormalized) variance of the pion norm 
cr 2 (n) = L 4 ■ Var(n) appears to be a suitable 'or- Figure 3. Phase diagrams for WA(a) and MA(b). 

der parameter'. 



Figure 2. a 2 (II) for qMA at (3 = 0.8. 



Fig.2 shows <7 2 (II) for qMA at (3 = 0.8. There 
is a clear signal at some k c . With increasing L 
this peak becomes even more sharp. 

For dynamical fermions the signal becomes less 
sharp (at least at small L) because the fermionic 
determinant suppresses small eigenvalues. How- 
ever also in this case a 2 (n) develops a pronounced 
maximum thus allowing to locate k c . 

In Fig. 3 we show the phase diagrams for both 
dWA and dMA. At the moment we have no clear 
interpretation of the phase in the upper left cor- 
ner in Fig. 3a. 



5. Conclusions 

For the standard compact Wilson action we 
observe a presumably 1st order PT which disap- 
pears after suppressing lattice artifacts. 

These lattice artifacts influence strongly the 
distribution of eigenvalues Xi of the fermionic 
matrix M. This influence is especially pro- 
nounced for the near-to-zero values of Ai. 

After suppression of artifacts a chiral transition 
(i.e., appearance of near-to-zero eigenvalues of 
M.) is left on a 'horizontal' line k = k c ((3). 

As a preliminary conclusion, we have no sign 
for a qualitative change of the behaviour along 
the line k = n c {j3). 
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